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FOREWORD  ..  V  . 

This  report  was  prepared  by  the  MRD  Division  of  the  General  American 
Transportation  Corporation  under  Contract  OCD-OS-63^176  for  the  Office  of 
Civil  Defense.  Mr.  Frank  C.  Allen  of  OCD's  Directorate  of  Research  was 
project  monitor.  This  report  'covers  the  period  of  June  1963  to  August  1964. 

The  contract  was  a  study  of  the.  habitability  of  typical  identified  shelter 
configurations  with  consideration  given  to  various  space-  allotments,  ventilating 
rates,  environmental  criteria, land  structural  particulars.'  All  of, these  parum- 
eters  were  studied  over  a  wide  range  of  variation  with  the  primary  emphasis  on 
aboveground  shelter  spaces. 


Analytical  studies  of  the  environment  in  shelters  are  continuing  under 
a  subcontract  with  Stanford  Research  Institute,  and  a  comprehensive  report 
will  be  prepared  at  the  conclusion  of  planned  work. 
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ABSTRACT 


Transient  and  steady  analysis  are  used  to  determine  the  psychrometric 
conditions  that  develop  in  large  aboveground  fallout  shelters  ventilated  with 
unconditioned  ambient  air.  Thebe  analyses  consider  the  shelter  size.,  geometry 
and  construction,  the  psychrometric  condition  of  the  ambient  weather,  and  the 
various  metabolic  and  nonmetabolic  heat  loads  to  the  shelter  air'.  The  results 
of  this  study  indicate  that  during  the  hot  s.ummer  weather,  only  a  small  fraction 
of  the  total  energy  input. to  the  shelter  is  lost  through  the  shelter  boundary  3ur 
faces.  Thus,  the  ventilation  requirements  for  large  aboveground  shelters  can 
be  obtained  by  the  use  of  an  analysis  which  neglects  the  heat  loss  through  the 
shelter  boundaries.  This  means  that  aboveground  shelter  ventilation  systems 
should  be  designed  to  remove  the  entire  thermal  load  generated  within  the  shelter 

INTRODUCTION 

The  fallout  shelter  must  protect  the'  shelter  occupants  from  radioactive 
fallout  and  from  the  detrimental  effects  of  excessive  carbon  dioxide  levels, 
insufficient  amount  of  oxygen,  and  excessive  shelter  temperatures  and  humidities. 
The  shelter  ventilation  system's  equipment  insures  that  these  additional  safe¬ 
guards  are  maintained.  For  example,  the  carbon  dioxide  level  in  a  shelter  can 
be  kept  at  an  acceptable  value  by  ventilating  the  shelter  with  3  cfm  per  occu¬ 
pant.  But,  as  much  as  60  efm  of  outside  air  could  be  required  to  keep  the 
shelter  temperature  and  humidity  down  to  tolerable  limits  during  the  hot  summer 
months.  In  some  situations,  ventilating  with  external  ambient  air  will  not 


sufficiently  control  the  shelter  effective  temperature  and  air  conditioning 
equipment  will  he  required. 

To  properly  select  the  ventilation  equipment,  the  shelter  designer  must 
accurately  know  the  .ventilation  rate  required.  The  major  parameters  in  this 
determination  are  the  number  of  shelter  occupants,  the  shelter's  geographical 
location,  and  the  shelter’s  size  and  construction.  The  number, of  occupants 
establishes  the  level  of  latent  and  sensible  metabolic  heat  that  is  produced 
in  the  shelter.  The  geographical  location  determined  that  statistical  nroba- 
bability  and  frequency  of  occurrence  of  the  psychrometric  condition  , of  the  ambient 
weather  from  which  the  weather  design  criteria  can  be  obtained.  A  knowledge 
of  the  shelter's  physical  size  and  construction  features  defines  the  heat 
transfer  coefficients  which. are  required  to  compute  the  heat  transfer  to  and  "  ■ 
from  the  shelter.  When  these  factors  and  the  required  shelter  effective  temper¬ 
ature,  are  known  and  incorporated  into  a  prediction  technique  for  the  ventilation 
requirements  of  the  shelter,  the  designer  has  a  basis  upon  which  he  can  select 
the  ventilation  equipment.  •  ... 

The  MED  .Division  has  been  engaged  in  the  development  of  such  a  prediction 
technique  under  the  Office  of  Civil  Defense  Contract  ODD -OS -63 -176,  Subtask 
1215A*.  The  eventual 'goal  of  this  program  is  to  formulate  a  simplified  procedure 
for  predicting  the  ventilation  requirements  of  shelters.  This  paper  presents 
the  results  of  several  studies  which  are  preliminary  steps  .in  the  development 
of  this  simplified  procedure.  These  studies  are  based  upon 

1.  the  metabolic  energy  output  of  approximately  400  Btu/hr- occupant 
given  by  the  sensible  and  latent  heat  expression,  see  Reference  1, 

•'(Mr.  F,  V,  Alien  of  O.  J' a  Directorate  of  Research  was  project,  monitor. 


2,  the  shelter  habitability  criteria  represented  by  the  ASHRAE  effective 
temperature  (an  empiracally  devised  index  of  the  various  phychrometric 
conditions  that  produce  similar  comfort  levels),  see  Appendix  A,  and 

3.  the  derived  mathematical  shelter  model  which  is  applicable  to  large 

.  shelters  that  are  either  aboveground  or  belowground,  do  or  do  not  have 
boundary  surface  heat  losses,  and  are  in  a  transient  or  steady-state 
condition  of  mass  and  energy  transfer,  see  Appendix  B  and  Reference  1. 

The  mathematical  model  for  the  shelter  is  the  most  comprehensive  analytical 
procedure  that  could  be  developed  without  unduly  complicating  the  computational 
procedure.  When  used  as  a  transient  analysis,  the  model  permits  the  psychrometric 
State  of  the  shelter  air  to  be  computed  as  a  function  of  time.  The  analysis 
considers 

1.  time- varying  inlet  air  conditions, 

2.  time  varying  energy  Inputs  to  the  shelter  air  from  equipment  and  lights, 

3.  time  varying  solar  loads  which  have  been  transmitted  through  windows 
into  the  shelter  (i.e.,  time  varying  values  of  transmitted  3olar 
radiation  as  obtained  from  the  ASHRAE  Guide  and  Data  Book) , 

4.  metabolic  latent  and  sensible  energy  loads  based  upon  the  instantaneous 
psychrometric  state  of  the  shelter,  see  Reference  1,  and 

5.  shelter  boundary  surface  heat  losses  (or  gains)  based  upon  a  one-dimensional 
heat  transfer  analysis  neglecting  corner  effects. 

The  shelter  model  is  based  upon  the  assumption  that 

1.  the  air  within  the  shelter  is  so  completely  mixed  that  all  of  the 
shelter  air  .is  at  one  psychrometric  condition, 

2.  the  convective  heat  transfer  coefficient  for  each  external  boundary 
surface  of  the  shelter  is  not  a  function  of  temperature  and  has  a 
constant  value  over  the  surface, 
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3.  the  walls  and  floors  which  are  internal  to  the  shelter  volume  are  at 
the  dry-bulb  temperature  of  the  shelter  air, 

4.  the  radiative  energy  interchanges  within  the  shelter  can  be  neglected, 

5*  the  solar  direct  and  indirect  radiative  input  to  the  shelter  through 

windows  and  the  equipment  and  lighting  heat  loads  in  the  shelter  can 
be  grouped  together  as  one  time  dependent  load  factor  which  is  termed 
the  thermal  load, 

6.  the  ventilating  air  exhausted  from  the  shelter  is  at  the  psychrometric 
condition  of  the  shelter  air, 

7.  the  effects  of  condensation  on  the  walls,  floors,  and  ceilings  of  the 
shelter  can  be  neglected,  ' 

8.  the  thermal-physical  properties  of  the  structural  materials  are  not 
temperature  dependent,  and 

9.  the  solar  radiation  absorption  on  opaque  shelter  bondary  surfaces 
can  be  neglected. 

Thus  the  shelter  is  idealized  as  an  enclosed  volume  in  which  sensible 
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such  a  model  are  the  conservation  of  energy  and  the  conservation  of  the  masses 
of  dry  air  and  water  vapor. 

The  primary  aim  of  this  paper  is  to  determine  the  most  simple  analytical  ■ 
model  which  can  predict  the  psychrometric  conditions  which  develop  in  above- 
ground  shelters.  To  evaluate  this  model,  actual  shelter  test  data  are  compared 
to  computed  data;  and  variation  of  parameter  approaches  are  applied  to  several 
types  of  shelters  to  generalize  the  results.  The  analytical  shelter  models 
studied  in  the  paper  are 

1.  aboveground  adiabatic  boundary  shelter  model  (no  heat  transfer 
through  the  shelter  boundary  surfaces)  with  time -dependent  parameters 
(l.e.,  the  state  of  the  inlet  air  and  the  thermal  loads  added  to  the 
shelter  air), 

2.  aboveground  nonadiabatic  boundary  shelter  model  with  time -dependent 
parameters , 

3.  aboveground  adiabatic  boundary  shelter  model  with  steady-state' 
(constant  with  time)  parameters, 

4.  aboveground  nonadiabatic  boundary  shelter  model  with  steady-state 

parameters.  ■  '? 

RESULTS  OF  SHELTER  ANALYSIS 
Transient  Analyses  of  Shelters 

The  time-history  of  the  temperature  and  humidity  within  a  shelter  is  pre¬ 
dicted  by  the  transient  analysis  which  considers  the  time  varying  mass  and 
energy  balance  about  the  shelter.  The  validity  of  this  analysis  is  established 
by  comparison  of  several  sets  of  analytical  computations  with  data  from 
actual  shelter  tests.  The  shelter  tests  chosen  for  the  comparison  are  the 
MRD  Wilmington,  North  Carolina  test  #7  (Reference  2),  the  University  of  Florida 


Central  Stores  Building  test  Phat.  .  IV  (Reference  3),  and  the  MRD  Houston,  Texas 
test  II  (Reference  4).  The  Wilmington  shelter  is  a  210-man  aboveground  shelter, 
the  Florida  shelter  is  a  250-man  basement  shelter,  and  the  Houston  shelter  is 
a  290-man  basement  shelter.  Using  the  observed  inlet  air  data  and  the  shelter 
dimensions  and  construction  demils,  the  transient  analysis  computed  the  dry- 
bulb  and  effective  temperatures  for  the  shelter  air  as.  a  function  fj.f  time,  see 
Figures  2,  3,  and  .4.  The  instantaneous  calculated  values  are  generally  within 
2°F  of  the  experimental  data  for  the  shelter  dry-bulb  temperature  and  within 
1°F  for  the  shelter  effective  temperature. 

When  the  shelter  is  assumed  to  have  adiabatic  boundaries  (no  heat  transfer 
through  the  shelter  boundary  surfaces),  the  analytical  and  experimental  results 
agree  on  the  24  hr  average  to  within  2°F  for  the  dry-bulb  temperature  and  within 
1°F  for  the  effective  temperature.  The  adiabatic  boundary  results  are  consis¬ 
tently  at  or  greater  than  the  experimental  values ,  This  means  that  all  of  the 
tested  shelters  lost  energy  through  their  boundary  surfaces  and  that  the 
transient  analysis  siightly  overestimates  the  amount  of  this  energy  loss. 
Steady-State.  Analysis  of  Shelters  with  Boundary  Heat  Loss 

The  time -average  psychrometric  condition  of  the  shelter  can  be  determined 
by  the  steady-state  analysis  which  considers  the  time-average  values  of  the 
psychrometric  condition  of  the  inlet  air  and  the  heat  loads  generated  within 
the  shelter.  The  heat  loss  from  the  boundary  surfaces  of  an  aboveground  shelter 
is  determined  by  the  temperature  difference  that  exists  between  the  shelter  air 
and  the  ambient  air  external  to  the  shelter  and  by  the  heat  loss  coefficient,  UA. 
The  UA  value  for  a  shelter  is  determined  by  the  size,  geometry,  and  composition 
of  its  walls,  floor,  and  ceiling;  and  its  value  is  determined  by 
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Time-Average  Temperatures  of  Shelter 
Environment* 

Dry-Bulb, Tdb  Effective, 

-  :  Experimental  Data  88.8°F  8U.5°F 

-  :  Adiabatic  Boundary  90.5°P  8b.6°F 

Computation 

Nonadiabatic  Boundary  88.3°F  84 , 5°F 

Computation 

*For  20.0  CFM/OCCUPANT  Ventilation  Rate,  83.9°F  Average  Inlet  Air  Dry-Bulb 
Temperature  and  80'’F  Average  Inlet  Air  Effective  Temperature. 

FIGURE  2  -  COMPARISON  OF  COMPUTATIONS  AND  EXPERIMENTAL 
DATA  FOR  WILMINGTON,  N,C„  TEST  7 
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Time-Average  Temperatures  of  Shelter 
Environment* 

Dry-Bulb,  T^B  Effective, 

-  :  Experimental  Data  88.6°F  83.2"F 

-  :  Adiabatic  Boundary  91.5°P  84.2°F 

Computations 

;  Nonadiabatic  Boundary  87.5nF  82.0^ 

Computations 

*For  13.9  CFM/OCCUPANT  Ventilation  Rate,  8l.5"F  Average  Soil  Temperature, 
85°F  Average  inlet  Air  Dry-Bulb  Temperature  and  y8.3°F  Average  Inlet  Air 
Effective  Temperature 

FIGURE  3  -  COMPARISON  OF  COMPUTATIONS  AND  EXPERIMENTAL 
DATA  FOR  CENTRAL  STORES  BUILDING  GAINSVILLE, 

FLORIDA 
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Time-Average  Temperatures  of  Shelter 
Environment* 


Dry-Bulb,  T 

- - :  Experimental  Data  88.7°F 

-  !  Adiabatic  Boundary  Computation  90,5°F 

Nonadiabatic  Boundary  87.7°F 

Computation 


Effective,  T_ 
'  E 

84.8°f 

86.1°P 

84.5°F 


*For  12.8  CFM/OCCUPANT  Ventilation  Rate,  87.9°F  Average  Soil  Temperature, 
83.0°F  Average  Inlet  Air  Dry-Bulb  Temperature  and  79-°°F  Average  Inlet  Air 
Effective  Temperature 

FIGURE  4  -  COMPARISON  OF  COMPUTATIONS  AND  EXPERIMENTAL 
DATA  FOR  HOUSTON,  TEXAS  TEST  II 
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where :  F 

m  -  number  of  boundary  surfaces 

Uj  =  overall  heat  transfer  coefficient  of  boundary  surface  J* 

A  =  surface  area  of  boundary  surface  i 

U 

p  =  number  of  shelter  occupants 
To  evaluate  the  accuracy  of  the  analysis ,  the  Wilmington  aboveground 
shelter  test  data  are  compared  to  the  results  computed  by  the  steady-state 
analysis  with  boundary  heat  loss .  The  comparison  shows  that  the  time-average 
shelter  dry-bulb  temperature  is  determined  to  about  1.5°F  and  the  time -average 
shelter  effective  temperature  is  determined  to  about  2°F,  see  Table  1, 

TABLE  1 


Ventilation  Rate,  cfm/occupant 
Average  Inlet  Dry-bulb  Temp,,  °F 
Average  Inlet  Effective  Temp.,  °P 
Average  Inlet  Relative  Humidity 

Average  Shelter  Dry-bulb  Temp , ,  °F 
Average  Shelter  Effective  Temp.,  °F 


Values  Predicted  by  Steady-State  Analysis  with  Boundary  Heat  Loss 


1.  ABOVEGROUND  SHELTER 

TEST** 

Test  Nos . 

5 

7 

8 

10A 

9.0 

20.0 

7.0 

13.0 

83.1 

83.9 

77.5 

90.4 

78.6 

80.7 

74,4 

80.1 

0.68 

0.78 

O.72 

0.4l 

90.0 

■  88.8 

'  89.8 

94.9 

85.7 

81).. 5 

86,1 

85.4 

Shelter  Dry-bulb  Temp.,  °F 
Shelter  Effective  Temp.,  °F 


*As  obtained  from  ASHRAE  Guide  and  Data  Book  -  Fundamentals  and  Equipment, 
American  Society  of  Heating,  Refrigeration,  and  Air  Conditioning  Engineers,  Inc,, 
1963,  pps.  391-te8. 

**For  210-man  aboveground  shelter  of  Reference  2. 


These  results  are  similar  to  the  agreement  obtained  with  the  transient  analysis 
If  the  steady-state  shelter  values  of  test  #7  are  compared  to  the  time-average 
shelter  values  of  the  transient  study  for  test  #7,  the  results  are  within  2°F 
of  each  other.  This  shows  that  the  steady-state  shelter  condition  approximates 
the  time-average  shelter  condition  determined  by  the  transient  analysis. 

The  steady-state  analysis  with  boundary  heat  losses  was  used  to  determine 
the  psychrometric  condition  of  the  shelter  as  a  function  of  the  ventilation  . 
rate  per  occupant,  the  equipment  and  lighting  load  per  occupant,  the  psychro¬ 
metric  conditions  of  the  inlet  air,  and  the  heat  loss  coefficient,  UA,  see 
Figure  5,  6,  and  7.  The  shelter  dry-bulb  temperature  variation  due  to  changes 
in  the  relative  humidity  of  the  inlet  ventilation  air  is  negligible  (less 
than  0.01°F).  This  is  not  true  of  the  effective  temperature  of  the  shelter 
and  therefore,  inlet  air  relative  humidities  of  15  and  80$  are  presented. 

The  main  effect  of  the  heat  loss  coefficient,  UA,  is  to  decrease  the 
shelter  dry-bulb  and  effective  temperatures  as  the  UA  value  increases.  But, 
the  temperature  reduction  decreases  as  the  flow  rate  is  increased,  see  Figure  8 
Generally,  less  than  20$  of  the  total  energy  input  to  the  shelter  is  lost 
through  the  shelter  boundaries;  however,  this  percentage  can  become  as  high  as 
30-50$  for  external  ambient  dry-bulb  temperature  below  70°F  and  flow  rates 
below  20  cf m/occupant. 

Steady-State  Analysis  of  Adiabatic  Boundary  Shelters 

The  adiabatic  boundary  shelter  model  differs  from  the  nonadiabatic 
boundary  shelter  model  in  that  it  neglects  any  heat  loss  through  the  shelter 
boundary  surfaces.  This  is  not  an  unrealistic  assumption,  because  the  heat- 
losses  of  the  aboveground  shelter  can  be  a  small  percentage  of  the  total  heat 
input  to  the  shelter.  In  the  belowground  shelter,  a  quasi -adiabatic  boundary 
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Shelter  Dry-Bulb  Tempereture 


105 


100 


90 


85 


80 


75 


Neglecting  solar  absorption  on  walls. 

Equipment  and  lighting  Load: 

- o  mu/at-occuPANT 

-  50  BOT/HR* OCCUPANT 

- 100  btu/hr-occupant 


In  All  These  Calculations 
Metabolio  Loada*400  BTO/HR-OCCUPANT 


10  20 


JL 

30 


kO  50  60 


Ventilation  Plow  Rate,  F]_,  CFM/OCCUPANT 


FIGURE  5  -  NONADIABATIC  BOUNDARY  SHELTER  DRY-BULB  TEMPERATURES  FOR  ALL  INLET 
AIR  RELATIVE  HUMIDITIES  AND  UA  =  10  BTU/HR- “F-OCCUPANT 


Inlet  Air  Dry-Bulb 


Inlet  Air  Dry-Bulb 
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FIGURE  7  -  NONADIABA'TIC  BOUNDARY  SHELTER  EFFECTIVE  TEMPERATURES  FOR  80f0 
INLET  AIR  RELATIVE  HUMIDITY  AND  UA  "  20  BTU/HR-°F-OCCUPANT 
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Dry- Bulb  Teapera-ture 


condition  is  reached  when  the  wall  temperature  approximates  the  shelter  average 
dry-bulb  temperature.  For  example,  consider  the  Houston  basement  shelter  tests. 
In  these  shelter  tests  the  time-average  soil  temperature  was  8t.9°F  whereaB 
the  time-average  shelter  dry-bulb  temperatures  varied  from  82  to  8f°F.  When 
the  Houston  shelter  is  considered  as  an  adiabatic  boundary  shelter  by  .the 
steady-state  analysis,  the  computed  data  agree  with  the  test  data  to  within  26F 
in  dry -bulb  temperature  and  1,5  °F  in  effective  temperature,  sec  Table  2. 


TABLE  2 

SUMMARY  OF 

HOUSTON  SHELTER 

'TEST* 

Test  Nos. 

Observed  Data 

I 

II 

III 

IV 

Ventilation  Rate,  cfm/occupant 

9.25 

12.8 

9.25 

18.5 

Average  Inlet  Dry-bulb., Temp. ,  °F 

82 

83 

77 

79 

Average  Inlet  Effective  Temp.,  °F 

78 

79 

76 

77 

Average  Inlet  Relative  Humidity 
Average  Shelter  Dry -bulb  Temp.,  °F 
Average  Shelter  Effective  Temp.,  °F 


O.77  0.75 

89.5  88.7 


Average  Shelter  Effective  Temp.,  °F  87  84.8 

Values  Predicted  by  Steady-State  Adiabatic  Boundary  Analysis 


0.79 

84 

82 


Shelter  Dry-bulb  Temp.,  °F 
Shelter  Effective  Temp.,  °F 


90.9  90.1 
87.4  85.9 


In  test  II,  the  steady-state  and  transient  results  agree  to  within  Q.5°F  when 
the  shelter  is  considered  to  have  adiabatic  boundary  walls.  This  again  confirms 
the  relationship  of  the  time-average  values  of  the  transient  analysis  with  the 
adiabatic  boundary  steady-state  analysis. 

The  shelter's  dry-bulb  and  effective  temperatures  were  computed  by  the 
adiabatic  boundary  steady-state  analysis  for  the  same  range  of  parameters  used 
with  the  nonadiabatic  boundary  analysis,  see  Figures  9>  10,  and  11.  When 
these  data  are  compared  with  the  data  for  the  nonadiabatic  boundary  shelter, 
the  inlet  air  conditions  that  produce  high  rates  of  energy  loss  are  found  to 


*For  290-man  belowground  basement  shelter  of  Reference  4. 
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ADIABATIC  BOUNDARY  SHELTER  EFFECTIVE  TEMPERATURES  FOR 
15$  INLET  AIR  RELATIVE  HUMIDITY 


Inlet  Air  Dry-Bulb  Temperature 


occur  when  energy  loss  is  the  least  needed  to  control  the  shelter  environment. 

For  example,  at  inlet  dry-bulb  temperatures  below  TO°F,  the  adiabatic  boundary 
shelter's  effective  temperature  never  exceeds  85°F  except  at  flow  rates  of  6  cfm/ 
occupant  and  below,  see ' Figures .10  and  11. 

If  the  effective  temperature  index  is  accepted  as  the  shelter  habi- 
ability  criteria,  a  psychrometric  chart  can  be  used  as  a  means  of  determining 
the  relationship  between  the  psychrometric  condition  of  the  inlet  air  and  the 
ventilation  flow  rate,  Bee  Reference  5-  For  instance,  if  an  average  effec¬ 
tive  temperature  of  8  5  0  F  is  chosen  as  the  habitability  limit,  the  loci  of 
the  average  psychrometric  conditions  of  the  inlet  air  that  produce  an  85 6F 
effective  temperature  in  the  shelter  can  be  plotted  for  various  ventilation 
rates,  see  Figure  12. 

Equipment  and  Lighting  Load  =  Zero 
Metabolic  Load  ®  400  Btu/hr -occupant 
Steady-State  Adiabatic  Boundary 
Shelter  Model 


FIGURE  1.2  -  VENTILATION  RATES  REQUIRED  TO  MAINTAIN  AN  AVERAGE  EFFECTIVE 

TEMPERATURE  OF  85°F 


Local  Ventilation  Rates 


In  order  to  determine  the  ventilation  rate  required  for  any  giyen  shelter, 
weather  design  criteria  must  be  selected.  Several  weather  data  studies  are 
currently  available  upon  which  the  weather  design  criteria  could  possibly  be 
based  (e.g.,  see  References  6,  7,  and  8),  but  the  weather  design  criteria  to 
be  selected  for  a  shelter  must  be  based  upon  a  detailed  study  of  the  relation¬ 
ship  between  the  psychrometric  state  of  the  ambient  weather  and  the  effect  that 
ventilation  air  from  this  ambient  has  on  the  shelter  environment  when  actual 
weather  data  are  considered.  This  relationship  will  be  influenced  by  the  Interval 
of  time  chosen  for  the  comparison.  Possibly  none  of  the  available  weather  data 
studies  will  be  applicable  and  a  new  set  of  criteria  will  have  to  be  established . 
l‘he  MRD  Division  is  presently  engaged  in  a  program  for  the  Office  of  Civil 
Defense  to  study  the  relationship  of  periods  of  actual  weather  on  shelter  envi¬ 
ronment.  In  addition,  more  information  must  be  gathered  concerning  the  effect 
of  environmental  conditions  on  the  human  body;  particularly,  when  these  con¬ 
ditions  are  changing  with  time  and  the  body  is  under  a  high  level  of  emotional 
stress.  However,  this  paper  has  shown  that  the  regional  ventilation  requirements 
for  shelters  can  be  determined  once  the  weather  design  criteria  and  habitability 
criteria  are  known, 

CONCLUSIONS 

The  transient  analysis  is  able  to  predict  an  aboveground  shelter's  instan¬ 
taneous  dry-bulb  temperature  to  within  2°F  and  the  shelter's  instantaneous 
effective  temperature  to  within  l"F,  The  time -average  values  of  a  shelter's 
dry-bulb  and  effective  temperature  are  predicted  by  a  steady-state  analysis  to 
within  2°F.  Either  of  these  analyses  can  be  used  to  determine  the  environ¬ 
mental  condition  of  a  shelter  with  a  reasonably  high  degree  of  accuracy.  The 
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steady-state  analysis  has  established  that  the  energy  loss  that  can  occur 
through  the  boundary  surfaces  of  an  aboveground  shelter  is  generally  less  than 
20$  of  the  total  thermal  energy  introduced  into  the  shelter  air  during  hot 
summer  weather.  All  of  these  results  indicate  that  the .mechanism  of  heat  loss 
through  the  boundary  surfaces  of  an  aboveground  shelter  cannot  be  depended 
upon  to  remove  energy  from  the  shelter  during  hot  weather.  At  most,  boundary 
surface  heat  loss  should  be  regarded  as  a  possible  safety-factor  in  ventilation 
system  design.  It  is  therefore  recommended  that  the  ventilation  systems  for 
aboveground  shelters  be  designed  to  remove  the  entire  thermal  load  generated 
within  the  shelter .  , 

The  flow  rate  predicted  by  the  analysis  of  a  shelter  without  boundary 
surface  heat  loss  is  the  maximum  flow  rate  that  can  be  required  (assuming  the 
solar  radiation  absroption  effects  on  opaque  shelter  boundary  surfaces  are 
negligible*) 5  and  thus  provides  a  means  for  establishing  an  upper  limit  on  the 
sise  of  the  ventilation  equipment  for  an  aboveground  shelter.  The  reasonable 
agreement  between  the  calculations  for  a  shelter  without  boundary  surface  heat 
loss  and  the  shelter  test  data  insures  that  the  equipment  size  based  upon  this 
ventilation  rate  is  not  overly  conservative. 

The  paper  has  shown  that  with  metabolic  head  load  data  and  weather  design 
criteria  the  ventilation  requirements  for  aboveground  shelters  (e.g,,  the 
shelters  surveyed  in  the  National  Fallout  Shelter  Survey)  can  be  determined 
analytically.  The  reliability  of  these  predictions  is  primarily  dependent  upon 
the  reliability  of  the  metabolic  and  weather  design  data  used  in  the  calculation 
procedure . 

*This  assumption  is  presently  under  detailed  evaluation  by  the  MRD  Division. 
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DERIVATION  OF  EFFECTIVE  TEMPERATURE  EQUATION 


Many  sets  of  psyehrometric  conditions  provide  the  human  body  with  similar 
levels  of  physiological  comfort.  The  effective  temperature  scale  has  been 
established  as  an  empirically  derived  index  of  the  various  psyehrometric  condi¬ 
tions  which  produce  similar  eomfort  levels.  The  effective  temperature  is  a 
function  of  the  dry-bulb  temperature,  wet-bulb  temperature,  and  velocity  of  the 
air  in  which  a  person  resides.  Generally,  these  data  have  been  presented  in 
nomographic  form  or  by  approximating  equations.  For  the  purposes  of  this  study, 
the  nomograms  were  inconvenient  to  use  since  they  were  not  in  equation  form, 
and  the  existent  equations  were  not  of  sufficient  accuracy.  As  a  result,  an 
accurate  mathematical  expression  for  effective  temperature  was  derived. 

The  derivation  was  based  upon  the  ASHRAE  nomogram*  for  effective  temperature, 
see  Figure  Al.  Because  the  ventilating  air  velocities  in  shelters  is  generally 
low,  only  the  20  ft/min  flow  velocity  curve  of  the  nomogram  was  considered. 

This  eliminated  the  a.ir  velocity  as  a  parameter  in  the  determination  of  effec¬ 
tive  temperature.  A  schematic  of  the  curve  AA  which  represents  the  20  ft/min 
flow  velocity  in  the  nomogram  is  shown  in  Figure  A2.  Curve  AA  is  linearized 
by  a  straight  .line  which  intersects  the  wet-bulb  (T,m)  and  dry-bulb  ('I')  temper- 
ature  scales  at  the  points  M  and  N  respectively.  By  geometric  similarity 


P-Q  M-Q 

N-0  '  M-0 


(Al) 


and 


P-Q  R-Q 

S-T  R-T 


(A2) 


*AGHRAE  Guide  and  Data  Book  -  Fundamentals  and  Equipment,  American  Society  of 
Heating,  Refrigeration,  and  Air  Conditioning  Engineers,  Inc.,  190.3,  p.  Ill 
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Substituting  the  above  values  into  Equation  (A3)  yields 


EFF 


107 (IDB 

62.3  +  (T 


•  V  * 62-3  <W 


DB 


(A4) 


This  relationship  is  limited  to  low  air  velocities  and  is  restricted  to  the 


temperature  range  of 

'.5-P  <  TM  <  110-F 


and 

30 °F  <  Twb  <  100 °F 

The  temperature  range  restrictions  are  imposed  due  to  the  fact  that  the  T-,™ 

ii).1  i 

curve  of  the  nomogram  cannot  be  considered  a  straight  line  beyond  these  temper¬ 
ature  limits.  Table  A1  shows  that  Equation  (A4)  has  an  average  error  of  less 


than  one  percent. 

TABLE  A1  ERROR  IN  CALCULATED  ERRECTIVE  TEMPERATURES  BY  EQUATION  (A4) 


T  .  °F 
DB’ 

V  °F 

Teff’  Cai'd  by 
Equation  (A4) 

Teff’  read  from 
Figure  (Al) 

Percent  Error 

110.0 

98.0 

99.4 

99-0 

0.45 

110.0 

90.0 

94.2 

94.0 

0.21 

110.0 

30.0 

73.5 

74.0 

'  0.68 

90.0 

80.0 

03.9 

83.9 

0 

90.0 

60.0 

75-5 

75.6 

0.13 

90.0 

4o.o 

70.0 

70.0 

0 

80.0 

60.0 

71.5 

71.5 

0 

80.0 

4o.o 

66.3 

66.3 

0 

80.0 

30.0 

64.7 

64.3 

0.62 

70.0 

60.0 

66.6 

66. 5 

0.15 

70.0 

50.0 

63.9 

64.0 

0.15 

1+5.0 

30.0 

45.1 

45.0 

0.22 

The  effective  temperature  values  are  correlated  to  the  relative  strain 
indicies  with  a  0.5  relative  strain  index  being  essentially  an  85"F 
effective  temperature,  see  Figure  A3-  Almost  all  people  are  comfortable  at 
zero  relative  strain  index  and  physical  failure  is  rapid  and  severe  at  1.0 
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relative  strain  index.*  This  comparison  illustrates  the  convertabi lity  of  the 
effective  temperature  index  to  other  indices  of  physiological  comfort. 


*11.  11.  Lee  and  A.  llensehel  "Evaluation  of  Thermal  Environment  in  Shelters",  U.S. 
Department  of  Health,  Education,  and  Welfare,  TR-8,  August,  1963 
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APPENDIX  B 


ANALYTICAL  MODEL 


B.l  DERIVATION  OF  ANALYTICAL  SHELTER  MODEL 


The  analytical  model  for  a  shelter  is  formulated  with  the  philosophy  that 
whether  a  shelter  is  above  or  below  ground  and  whether  it  is  single  story  or 
multi-story,  the  shelter  can  be  represented  by  a  single  mathematical  model. 
Furthermore,  the  model  is  developed  with  the  idea  that  it  should  be  the  simplest 
model  that  can  be  devised  and  still  adequately  explain  the  phenomena.  In  accord 
with  this  philosophy,  the  assumptions  are  made  that 

1.  the  air  within  the  shelter  is  so  completely  mixed  that  all 
of  the  shelter  air  is  at  one  psychrometric  condition, 

2.  the  convective  heat  transfer  coefficient  for  each  external 
boundary  surface  of  the  shelter  is  not  a  function  of  tempera¬ 
ture  and  has  a  constant  value  over  the  surface, 

3.  the  walls  and  floors  which  are  internal  to  the  shelter  volume 
are  at  the  dry-bulb  temperature  of  the  shelter  air, 

4.  the  radiative  energy  interchanges  within  the  shelter  can  be 
neglected, 

5.  the  solar  direct  and.  indirect  radiative  input  to  the  shelter 
through  windows  and  the  equipment  and  lighting  heat  loads  in  the 
shelter  can  be  grouped  together  as  one  time  dependent  load  factor 
which  is  termed  the  thermal  load, 

6.  the  ventilating  air  exhausted  from  the  shelter  is  at  the 
psychrometric  condition  of  the  shelter  air, 

7.  the  effects  of  condensation  on  the  walls,  floors,  and  ceilings 
of  the  shelter  can  be  neglected,  and 

8.  the  thermal -physical  properties  of  the  structural  materials 
are  nob  temperature  dependent. 
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The  model  of  the  shelter  is  established  by  these  assumptions.  Thus  the  shelter 
is  idealized  as  an  enclosed  volume  in  which  sensible  heat,  latent  heat,  solar 
heat,  equipment  heat,  and  ventilating  air  are  introduced' and  from  which  air  is 
exhausted  and  energy  is  lost.  The  governing  principles  for  such  a  model  are 
the  conservation  of  energy  and.  the  conservation  of  the  masses  of  dry  air  and 
water  vapor . 


The  conservation  of  energy  requires  that 


the  rate  of  change  of  enthalpy  the  rate  of  c 

in  the  ventilating  air  +  in  the  shelter  air 


J 


the  rate  of  metabolic  heat 
input  by  the  shelter  occupants 


j  ["the  rat 

iS  +  I  load  to 

and  integrating  Equation  (Bl)  over  the  time  increment  At 


the  rate  of  change  of  enthalpy-] 


ate  of  thermal”] 
to  the  shelter 


' 


the  rate  of  energy  transfer  across" 
the  shelter  boundary  surfaces 


(Bl) 


(“l  -  V  •  V>  *  S.  ♦  «T  -  %  -  0  <K> 

whore : 

-  enthalpy  of  ventilating  air  entering  the  shelter,  Btu 

Hg  =  enthalpy  of  ventilating  air  leaving  the  shelter,  Btu 

H„  .  =  enthalpy  of  shelter  air  at  beginning  of  time  increim  Ai  ,  if  o 

b,  1 

H0  „  =  enthalpy  of  shelter  air  at  end  of  time  increment  A".'  ,  btu 
b,  d 

-  energy  due  to  metabolism  of  shelter  occiipants,  Btu 

~  energy  input  due  to  thermal  load  (see  assumption  #5.>  p.  B2),  Btu 

-  energy  loss  or  gain  through  shelter  boundary  surfaces,  Blu 
a 


The  conservation  of  the  mass  of  water  vapor  necessitates  that 


t 

\ 


the  rate  of  change  of  water-] 


vapor  in  the  ventilating  nirj  +  j^vapor  1 


Ptbe  rate  of  chang 


in  the  shelt 


;e  of  water 
Lter  air 


the  rate  of  water  vapor 

introduced  by  the  shelter  oceupantsj  =  0 


s 


and  integrating  Equation  (B3)  over  'the  time  increment  At 


(Vl  '  MV,2)  +  (\s,l  '  ^,S,2)  +  *V/] 


(B3) 


(bU) 


where : 


My  j  -•  mass  of  water  vapor  in  the  entering  ventilating  air,  lb 

My  g  =  mass  of  water  vapor  in  the  exhausting  ventilating  air,  lb 

My  g  ^  =  mass  of  water  vapor  in  the  shelter  at  the  beginning  of 

’  ’  the  time  increment  kr  ,  lb 

My  p  -  mass  of  water  vapor  in  the  shelter  at  the  end  of  time 

’  ’  increment  A  T  ,  lb 

My  M  «  mass  of  water  vapor  introduced  by  the  shelter  occupants,  lb 
The  conservation  of  the  mass  of  dry  air  demands  that 


the  ra 
air  in 


rate  of  change  of  dry 
the  ventilating  air 


y  the  rate  of  change  of  dry 

ir  +  air  in  the  shelter  air  •- 


and  integrating  over  the  time  increment  A r 


(M  .  -  M  )  +  (M  -  M  _  c 

a,l  a, 2  a,S,l  a,S,2 


)  =  0 


0 

(B5) 

(b6) 


where: 


a,  1 
^a,  2 


M 


a,  S,  1 


a ,  o ,  2 


mass  of  dry  air  in  the  entering  shelter  air,  lb 

mass  of  dry  air  in  the  exhausting  shelter  air,  lb 

mass  of  dry  air  in  the  shelter  at  the  beginning  of  the 
time  increment,  lb 

mass  of  dry  air  in  the  shelter  at  the  end  of  the  time 
increment,  lb 


Denoting 


AQ  =  li±  +  HSjl  +  Qn  +  -  Qb  (B7) 

Vo  =  Vl  +  V  S,  1  +  V  M  (b8) 

Ma,0  =  Ma,l  +  Ma,S,l  (B9) 

and  substituting  Equation  (B7)  through  (B9)  into  Equations  (B2),  (b4)  and 
(B6)  gives. 

,  AQ  =  Hg^g  +  Hg  (BIO) 

^,0  "  \2  +  \S;2  (Bll) 

Ma,0  =  Ma,2  +  Ma,S,2  (B12) 

where i 

HS,2  =  Ma,S,2ha,S,2  +  Vs,2\s,2  (B13) 

,H2  =  Ma,2ha,2  +  V2hV,2  (Bl4) 

From  the  assumption  that  the  exhaust  air  has  the  same  psychrometric  condition 
as  the  shelter  air 

ha,2  =  ha,S,2  (B15) 

and 

hV,2  \s,2  (Bl6) 

With  Equations  (Bll)  to  (B1.6)  substituted  into  Equation  (BIO) 

=  Ma,0ha,S,2  +  VohV,S,2  (B17) 

The  specific  enthalpies  of  dry  air  and  water  vapor  are  given  by 

h  =  0.24  (Tnn  +  459-69)  (B18) 

3,  19  Jj 

and 

h  -  1061  +  0.444  (T„J  (B19) 

v  .DB 

for  32°F  <  Tdb  <  150°F 

b:> 


where : 


=  dry-bulb  temperature,  °F 

Combining  Equations  (Bl?)  to  (B19)  and  solving  for  the  dry-bulb  temperature 
results  in 


Vo 

DB, 2  0.444  My  0  +  0.24  M&  ~ 


(B20) 


where : 


dry-bulb  temperature  of  shelter  air  at  end  of  time  increment 

At  ,  y 

The  following  relationships  are  evident: 


.  Ma,2  = 

6°P2pa,2AT 

(B21) 

^V,  2 

60  *2  p  V,  2  AT 

(B22) 

Ma,S,2 

V  Pa,S,2 

(B23) 

and 

\s,2  = 

V  PV,S,2 

(B24) 

where : 


V  =  shelter  volume,  ft 
At  =  time  increment,  hr 

O 

0  =  density  of  dry  air  leaving  shelter,  lb/ft J 

Qi)  cL 

,Oy  g=  density  of  water  vapor  leaving  shelter,  lb/ft3 


p  0  -  density  of  dry  air  in  shelter  at  end  of  time  increment, 

a’°’2  .lb/ft3 


P 


V,S,2 


t 


2 


density  of  water  vapor  in  shelter  at  end  of  time  increment 
lb/ft3 


volumetric  flow  rate  of  exhaust  air, 


ft3/min 


but  by  assumption 


/a,  2 

pa,  S ,  2 

(B25) 

and 

P  V,2 

'  PV  .  0  ,  r> 

(B26) 

Bo 

Substituting  Equation  (B2l)  through  (B26)  into  Equation  (Bll)  and  (B12)  yields 

\,0  "  PV,2  (6°  F2Ar  +  V)  (B2Y) 


and 


".,0  ■  <“  V  *V> 

When  Equations  (B27)  and  (B28)  are  solved  for  Eg,  they  result  in 

M  -  p  V 
a,  0  'a,  2 


5oSr7a>2 

ny  .o  ~pv,2v 


6oat  p. 


(B28) 

(B29) 

(B30) 


V,2 


The  densities  of  the  dry  air  and  water  vapor  are  assumed  to  obey  the  perfect 
gas  law.  Thus; 

P 


V 


and 


where : 


v  ~  Rv(Tdb  +  459.89) 


p  -  p 
B  V 

P  a  Ra(TDB  +  1+59.69) 


(■B31), 


(B32) 


PB  =  barometric  pressure,  lb/ft 

2 

Py  =  partial  pressure  of  water  vapor,  lb/ft 

Ry  =  gas  constant  for  water  vapor,  ft-lb/lbmole-°R 

R  =  gas  constant  for  dry  air.  ft-lb/lb  ,  -°R 

a  '  mole 

The  partial  pressure  of  water  vapor  is  expressed  by 


r  P„ 


(B33) 


where : 


relative  humidity 


P_  ■=  saturation  pressure  of  water  vapor,  lb/l'te 


Pg  --  '3 .  .1 32e°  ’ (  tbb  ) 


( B3ji : 


fur  32“f<  Tnp 


1'30°E 


BY 


which  is  an  analytical  expression  for  the  tabulated  values  of  saturated  vapor 
pressures  as  a  function  of  dry-bulb  temperature  (Ref.  8  ). 

From  Equation  (B30)  through  (B33)  with 


and 


»«  *  53-35  n-lb/1WR 


”v  -  85.71 


the  following  expressions  can  be  obtained: 

PV 


P„  a  r  [5.132e°‘329^TDB^] 


and 


P, 


l 


a  (,tdb  +  d>y.oy;  •  B 

Equating  Equation  (B27)  and.  (B28)  gives 


M. 


a,  0 


“v.O 


(B35) 


[  0.05987  r  e°'03i!9<V]  (B36) 

[O.Ol8?44P  -  0.0962  r  e° • °329(TDB) ]  (537) 


pa,2  PV,2 

Evaluating  Equations  (B36)  and  (B37)  at  the  shelter  dry-bulb  temperature, 
[tdb  gJ,and  substituting  the  results  into  Equation  (B38)  gives  the  partial 
pressure  of  the  water  vapor  [Py].as 


(B3B) 


(B39) 


v,o 


And  from  Equation  (B33)>  the  relative  humidity  of  the  shelter  at  the  end  of 


the  time  increment  is  given  by 


r  - 


(B4o) 


Substituting  Equations  (B34)  and  (B39)  into  Equation  (b4o)  results  in 


Bo 


V 


(B4l) 


5-132e  (0.0329(Tdb)>  [l  *  0.6221,  i°l 

1  “v,oJ 


With  these  equations,  the  psychrometric  conditions  within  the  shelter  can  be 
determined  as  a  function  of  time  through  the  use  of  the  following  procedure. 

1.  Determine  the  psychrometric  conditions  of  the  air 
introduced  into  the  shelter  and  the  air  within  the  • 

,  shelter  along  with  the  heat  inputs  and  losses  of 

the  shelter  for  the  time  interval  At 

2.  Compute  the  quantities  AQ,  M,,  Q,  and  Ma  Q  from  Equations 
(B7)>  (b8),  and  (E9)  respectiv4ly.  ’ 

3.  Assume  that  these  quantities  do  not  change  over  the  time 
interval  AT,  and  compute  the  dry-bulb  temperature  of  the 
shelter. at  the  end  of  the  time  interval  from  Equation  (B20) . 

4.  Calculate  the  relative  humidity  of  the  shelter  air  at  the 

end  of  the  time  interval  from  Equation  (B4l) .  \  . 

5.  Set  the  psychrometric  conditions  of  the  shelter  at  the  end  of 

the  time  increment  equal  to  those  in  the  shelter  for  the  beginning 
of  the  next  increment,  and  repeat  the  entire  procedure. 

By  the  continuous  application  of  this  procedure,  the  shelter's  psychrometric 

condition  can  be  obtained  as  a  function  of  time  for  any  time  period.  The 

computation  method  that  has  been  developed  constitutes  a  transient  analysis 

of  the  shelter  environment  under  the  influence  of  time  varying  parameters. 


The  shelter's  psychrometric  condition  is  described  in  the  analysis  by  its 
dry-bulb  temperature  and  relative  humidity.  Instead,  the  dry-bulb  and  wet-bulb 


or  effective  temperatures  can  be  used.  The  wet-bubl  temperature  is  determined 


from  the  Carrier  equation. 


l'PB  "  Ps'^TDB  ~ 

2800.0-1.3  (twb) 


(b42) 


J.  H.  Carpenter,  "Fundamentals  of  Psychrome tries" , 
Technical  Note  T  200-20,  Syracuse,  New  York,  1962, 


Carrier  Corporation, 
p.  12. 


B9 


where: 


2 

Pg1  =  saturation  pressure  at  wet-bulb  temperature,  lb/ft 

T,_  =  wet-bulb  temperature,  °F 

Wi3 

with  P  1  evaluated  by  Equation  (B34) 

Pg'  =  5-132  e0,0329^WB)  (BU-3) 


and  obtained  from  Equation  (B39) •  The  results  is  a  transcendental  equation 
that  can  be  solved  for  the  wet-bulb  temperature,  Tyg.  When  the  value  of  the 
wet-bulb  temperature  is  known,  the  effective  temperature  can  be  determined 


from  equation  (A4) . 


107.5  <tdb  -  W  +  &'3  <V 


=2'3  *  <tdb  -  W 


As  mentioned  in  Appendix  A,  this  relationship  is  limited  to  low  air  velocities 
and  is  restricted  to  the  temperature  range  of 

45eF  <  T'  <  110 nF 

—  JJij 

and 

30  °F  <Tra<  100  °F 


The  heat  loss  or  gain  of  the  shelter  boundary,  Q^,  over  a  time  interval  At 
is  determined  by  the  temperature  of  the  inner  surface  of  the  boundary,  T^.  In 
the  computations,  the  value  of  T^  is  obtained  from  the  temperature  distribution 
that  existed  through  the  boundary  during  the  previous  time  increment.  This 
temperature  distribution  is  deduced  from  a  transient  analysis  of  the  energy 
transfer  in  the  boundary.  In  order  t.o  accomplish  this  analysis,  the  boundary 
is  divided  into  a  finite  number  of  slabs  of  thickness  AX,  except  for  the  inner 
and  outer  surfaces  which  are  made  into  slabs  of  thickness  A  X/2,  with  each  slab 
assumed  to  be  at  a  single  temperature.  This  temperature  is  assigned  to  the 
midpoint  of  each  slab,  except  for  the  innermost  and  outermost  slabs  which  are 


assigned  temperatures  at  the  external  surfaces  of  the  slab  (see  Figure  B.l ) . 
These  locations  of  slab  temperature  are  termed  nodal  points.  The  external 
nodal  points  of  the  boundary  transfer  energy  by ' convection  with  surrounding  air 
arid  by  conduction  with  the  next  internal  nodal  poipt.  All  of  the  rest  of.  the 
nodal  points  transfer  energy  by  conduction  with  the  two  adjacent  nodal  points. 
The  energy  balance  about  an-  internal  nodal  point  in  is 


"energy  conducted"" 

"energy  conducted  to 

[""energy  stored  in  ] 

to  point  m  from 

...  - 

point  (m+1 )  to 

-• . . 

[the  slab  of  point  ml 

point  (m-l) 

point  m 

(34*3) 

(T  -T  J 
v  m  m-1' 


(T  ^ -T  ) 
m+1  nr 


AX>C  (T  '-T  ) 
pm  m 


and  . the  temperature  at  the  nodal  point  is 


(B47) 

(b48) 


k  =  thermal  conductivity  of  composite  boundary  material, 
Btu/hr-ft-°F 

y  '  o 

p  -  density  of  composite  boundary  material,  lb/ft J 

C  -  specific  heat  of  composite  boundary  material,  Btu/lb- °F 

Ax  -  thickness  of  boundary  slab,  It 

T  -  temperature  of  nodal  point  in  at  beginning  of  time 
111  increment,  "F 

T  ^  =  temperature  of  nodal  point  m-.'l  at  beginning  of  time 
increment,  ‘‘F 


T  =  temperature  of  nodal  point  m+1  at  beginning  of  time 
m  increment,  °F 

T  1  =  temperature  of  nodal  point  m  at  end  of  time  increment,  °F 

The  energy  balance  about  the  inner  surface  of  the  boundary  is 


energy  convected 
from  shelter  air 
to  boundary 


energy  conducted  to 
nodal  point  2  from 
_nodal  point  1 


energy  stored  in 

the  slab  at 

nodal  point  1  (b49) 


'  (t2-t  ) 

hiA(TDB,2  *  Tl>  + 


f-  A^(Tl'  '  Tl) 


and  the  temperature  at  the  inner  boundary  surface  is 


Ti  =  ■  Wl»,2  +  B3T2  +  (1'V  B2B3)T1 


with  B^Bg+l)  <  1 


where : 


hi  Ax 

k 


B„  =  2B, 


h  =  heat  transfer  coefficient  on  the  inner  surface  of  the  . 

1  boundary,  Btu/hr-ft2-°F. 

T,,,,  -=  temperature  of  shelter  at  beginning  of  time  increment,  °F 

JUBj  d 

=  temperature  of  nodal  point  at  inner  surface  of  boundary  at 
beginning  of  time  increment,  °F 

T  =  temperature  of  nodal  point  2  at  beginning  of  time  increment, 

2  oF 

'  =  temperature  of  nodal  point  at  inner  surface  of  boundary  at 
end  of  time  increment,  °F 

The  energy  balance  about  the  outer  surface  of  the  boundary  n  is 


energy  conducted 
from  (n-l)  nodal 
point  to  outer 
surface  point  n 


na.1  ambient  air 


energy  convected  energy  absorbed  energy  stored 

from  outer  surface  from  solar  radi-  in  the  slab 
point  n  to  exter-  +  ation  in  the  slab  --  at  the  outer 
nal  ambient  air  at  the  outer  ^surface 

surface 


surface 


BJ.3 


or 


-T-~— -  -  ‘.‘'WWb 


AX 


and  the  temperature  at  the  outer  boundary  surface  is 

T  ' 


where : 


At 

B, 


Vn-l  *  V3\  *  ^'h'hhK  *  % 


With  B3(B4  +  l)  <  1 


B,  = 


h0A* 


(B54) 

(B55) 

(B56) 


h 

o 


n-1 


=  heat  transfer  coefficient  on  the  outer . surface  of  the 
boundary,  Btu/hr-ft-  -°F 

=•  intensity  of  solar  radiation  for  shelter  location  and 
given  calender  date  per  unit  area,  Btu/hr-ft^ 

=  temperature  of  nodal  point  n  at  outer  surface  of  boundary 
at  beginning  of  time  increment,  °F 

temperature  of  nodal  point  n-1  at  beginning  of  time 
increment ,  °F 


T.  "  =  temperature  of. external  ambient  air  at  beginning  of  time 
increment,  °F 

T  '  =  temperature  of  nodal  point  n  at  outer  surface  of  boundary  ■ 
at  end  of  time  increment,  °F 

f  =  absorptance  of  outer  slab  surface  for  solar  radiation 

The  boundary  surface  heat  loss  as  defined  by  Equations  (b47),  (B5l),  and 
(B55)  is  applicable  to  the  one  dimensional  heat  transfer  from  a  boundary  sur¬ 
face.  If  the  boundary  surfaces  of  a  shelter  are  large,  corner  effects  are 
negligible  and  each  boundary  surface  ce.n  be  considered  to  be  conducting  energy 


from  or  to  the  shelter  uni-directionally;  However,  the  problem  in  solving  for 
the  heat  energy  loss  or  gain  from  the  shelter  boundary  surfaces  is  dependent 
upon  the  designation  of  a  temperature  in  Equation  (B55)  and  the  properties 
of  the  boundary  surface,  i.e.,  density  ( p  ),  thermal  conductivity  (k),  specific 


heat  (C  ),  and  thickness  (  X)  with  Av  =  — — 
v  p"  '  '  A  n-1 


where  n  =  number  of  nodal  points. 


To  simplify  the  analysis,  the  shelter  boundary  surfaces  are  grouped  together 
to  form  one  slab  surface  exposed  to  the  temperature  of  the  shelter  on  one  side, 
and  the  temperature  on  the  other  side.  The  procedure  of  replacing  all  of 
the  separate  boundary  surfaces  by  one  slab  surface  is  based  upon  the  determination 
of  area  weighted  average  values  for  the  various  properties  of  the  boundary  sur¬ 
face.  The  property  values  of  each  boundary  surface  are  based  upon  thickness 
weighted  average  property  values  of  the  materials  that  compose  each  of  the 
boundary  surfaces.  For  example,  for  any  boundary  surface  property  X,  the  area 

weighted  average  value  of'  property  X  for  the  slab  surface  is  given  by 

r(j)  m  A 

X  -  [  2  2  {^L)  (XpXp)  ]  /  A 


P-1 


(B57) 


where  2 


A. 

J 


-  surface  area  of  boundary  surface  ,j 
m  =  number  of  boundary  surfaces 

Xp  =  property  X  of  the  material  p  in  the  boundary  surface  j 
X  p  =  thickness  of  the  material  p  in  the  boundary  surface  j 
r(,-))  =  number  of  materials  in  boundary  surface  j 

r(j) 


X 


J 


A  = 


2 

P“1 

m 

2 

j-1 


X 


so  that 


r(j) 

[  S 

P-1 

r(j) 

[  s 


A. 


2 

J=1 

m 

2 

j=l 


A  O  1/  i 


P  P 


(  Xpkp)]  /A 


(B58) 


(B59) 


B15 


m 


r(j)  m 

[  2  2 

p=l  j=l 


(-jI)  (  A  (  C  )  , )  ]  /  A 
X  •  p *•  pJ  J  J  ' 

J 


( BuO ) 


and 


[2  A  \  ,  ]/  A 

J--1  J  J  ■ 


(B6l) 


The  values  p  ,  k,  C  '>  and  X  determine  the  property  values  of  the  composite 
boundary  surface  slab . 


In  general^  the  exterior  sides  of  shelter  boundary  surfaces  are  exposed 

to  three  types  of  environments.  The  boundary  surface  can  be  exposed  to  the 

ambient  weather  (e.g.,  an  outside  wall),  to  the  soil  (e.g.,  an  underground 

wall,  the  floor),  or  to  a  space  interior  to  the  structure  but  exterior  to  the 

shelter  (e.g.,  a  first  floor  shelter  ceiling  exposed  to  a  second  floor  space 

of  the  building).  In  each  of  these  situations,  the  value  of  T.  for  fitch 

A 

boundary  surface  will  be  different.  In  general,  the  area  weighted  average 

value  of  T,  would  be  . 

A  » 


where: 


2  a;(t.). 

d  A'j 

T  ~  - - - 

A  m 

2  A. 


temperature  to  which  the  exterior  side  of  boundary 
surface  j  is  exposed 


(B62) 


Letting 


P, 


m 

2  A. 
3 


with  j  -  1,2, ... , 


jl 


(B63) 


Blh 


m 


then 


T 


A 


m 

2 

J«1 


(B64) 


With  the  fact  that  the  exterior  side  of  the  boundary  surface  can  be  exposed 
to  only  three  types  of  environments,  then 

*A  =  Wo  +  ?3(Vs  +  Wi  vith  1  =  ,(B65) 

where: ' 

Pq  =  percentage  of  shelter  boundary  surfaces  exposed  to 
the  ambient  weather  temperature  (T^)q 

p  =  percentage  of  shelter  boundary  surfaces  exposed  to 
solid  at  the  temperature  (T^)g 

=  percentage  of  shelter  boundary  surfaces  exposed  to 
each  of  the  interior  space  temperatures  (T^)^  for 
each  of  the  q  interior  spaces. 


However,  experience  and  experimental  test  results  have  shown  that  the 


temperature  values  of  the  soil  and  interior  areas  adjacent  to  the  shelter  can 

-11  •»  . 

be  closely  approximated  by 


and 


(Vs 

(Vi 


Ts  +  tse 


(Vo  +  tse 


(b66) 

(b67) 


where: 


Tc  -  well. water  temperature  at  location  of  interest  as  given  by  Collins.* 

Tgg  -  shelter  dry-bulb  temperature  at  time  for  which  boundary  surface 
heat  loss  or  gain  is  computed. 

Substituting  Equations  (b66)  and  (B67)  into  Equation  (B65)  gives 


^A  '=  Po'Vo  f  2  [TS  +  Vi]  h  2  ['Vo  +  TSe] 


(B68) 


-K 


W.  D.  Collins,  "Temperature  of  Water  Available  for  Industrial  Use  in  the 
United  States",  U.S.  Geological  Survey  Water  Supply  Paper  520-F,  .1925. 


ty-r 

®e 


I 


i 


or 


..vhere: 


T  =  f-p  +  — )  (t  )  + 

A  2'  ''ho 


(ps  +  Pj) 


TSE  +  2  PSTS 


(B69) 


Pj  =  S  P^>  total  percentage  of  shelter  boundary  surfaces 
,1=1  .  exposed  to  interior  spaces 

The  value  of  f ^  of  Equation  (B69)  is  used  as  the  value  of  in  Equation  (B55) . 
With  these  area  weighted  average  property  values  and  the  relationship  for 
the  value  of  the  boundary  surface  heat  loss  or  gain,  Qg,  can  be  computed. 

The  analysis  also  considers  the  solar  radiation  that  enters  the  shelter 
through  windows  and  the  solar  radiation  that  is  absorbed  on  the  exterior  surfaces 
of  the  shelter r  The  solar  energy  that  comes  through  the  windows  is  treated  as  a 
time  varying  thermal  load  in  the  term  of  Eq.  (B2).  The  magnitude  of  this  load 
is  determined  by  ■'  | 

<V«inta,~  *  ■%>!  i}  W°> 


where  ■, 


(Qg) window  =  enerSy  transmitted  into  shelter  at  time  of  consideration 
«  solar  intensity  on  window  j  at  time  of  consideration 


7 


J 


=  surface  area  of  window  j 

=  transmittance  of  window  j  for  solar  radiation 
The  solar  energy  that  is  absorbed  by  the  walls  is  neglected  as  an  assumption 
in  the  study  of  the  shelters  mentioned  in  the  main  body  of  this  report.  Therefore 
in  Eq.  (B5*0,  i/j  is  taken  to  be  equal  to  zero.  Generally,  this  may  be  handled  in 
the  following  way. 

The  solar  energy  that  is  absorbed  by  the  wall  is  defined  by  the  solar  intensity 
Qg  times  the  absorptance  of  the  exterior  surfaces  for  solar  radiation,  ^  [  see  Eq. 
(B54)  ] .  The  value  of  that  is  used  in  Eq.  (B5^-)  is  an  area-weighted  average 

value  of  the  solar  intensities  that  occur  upon  each  of  the  boundary  surfaces  of 


BJ.O 


the  shelter  exposed  to  sunlight.  Consequently,  this  solar  energy  term  is  a 
function  of  time  and  is  defined  as 

i  w/j 


Z  A. 


where: 


A.  =  area  of  the  bondary  surface  J,  ft 
J 

n  =  number  of  bondary  surfaces  exposed  to  solar  radiation 
at  time  for  which  §g  is  required 

(Q  )  =  intensity  of  solar  radiation  incident  upon  the  bondary  2 

s  **  surface,  J,  at  time  for  which  Q  is  required,  Btu/hr-ft 


The  value  of  (Q  )  is  defined  by 
8  J 


«.),  -  «, 


direct  ^  j  +  ^diffuse  ^J  +  ^reflective 


where : 


(Q,,  .  )  .  =  energy  received  by  surface  j  directly  from  sun, 

direct  J  Btu/hr/ft2 

^diffuse  =  energy  from  atmospheric  scattering  which 

J  reaches  surface  J,  from  all  directions  in 

sky,  Btu/hr-ft^ 

(^reflective  K  =  enersy  reflected  from  objects  facing  the 
J  surface  jy  that  reach  the  surface  J, 

Btu/hr-  fir' 

B .  2  COMPUTATIONAL  CRITERIA 


The  time  increment  used  in  a  computation  cannot  be  arbitrarily  chosen  with¬ 
out  introducing  instability  into  the  computations.  Practice  has  shown  that  if 
the  volume  of  the  ventilating  air  introduced  into  the  shelter  during  a  time 
increment  is  less  than  a  fixed  percentage  of  the  internal  shelter  volume,  the 
calculations  will  be  stable.  That  is, 


B19 


sfeg  Ti  nflfcji  r.sESi  ■rx.-TsxxiKSXt*  *s 


60  (F  )  (At)  <  C  V 


where: 


fe(F1) 


C  =  a  constant  coefficient  (=  0.1,  by  experience)  . 

At  =  time  increment,  hr 

F^  =  volumetric  flow  rate  of  entering  ventilation  air,  ftVmin 

■5 

V  =  shelter  volume,  ft 


When  a  shelter's  stability  conditions  are  computed,  the  above  stability 
criterion  is  the  relationship  that  requires  the  smallest  time  increment,  At, 
and  not  the  relationships  of  Equations  (B48),  (B|52),  or  (B5^)  *  Therefore,  the 
Bize  of  the  time  increment  is  determined  by  Equation  (B7^)* 

\  ■  ^  <' 

Two  other  computational  procedures  that  are  followed  a\ /the  mathematical 
model  of  the  shelter  are 

1)  the  input  data  are  linearly  interpolated  for  time  increments  that 

,,  are  smaller  than  those  for  which  the  data  are  given,  and 

2)  the  relative  humidity;  of  Ijhe  shelter  is  always  kept  less  than  or  '■ 
equal  to  unity.  If  the  computed  relative  humidity  is  greater  than 
unity,  the  latent  heat  input  to  the  shelter  is  reduced  until  the 
relative  humidity  is  equal  to  unity, 

B3  -  STEADY-STATE  SHELTER  ANALYSIS 

In  the  steady-state  analysis  of  the  shelter  environment,  none  of  the  para¬ 
meters  vary  with  time  including  the  psychrometric  condition  of  the  shelter. 

This  reduces  the  complexity  of  the  computation  program.  The  transient  calcu¬ 
lation  procedure  can  be  used  to  obtain  steady- state  environmental  results  if  all 
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of  the  parameters  are  made  constant  with  time'  T!  process  is  an  iteration 
procedure  which  converges  to  the  steady- 8tat«  Values.  The  steady- s-ta^  caicu. 
latlon  Is  otherwise  identical  to  Use  tran#i@nt  shelter  model. 

Under  steady-state  conditions  the  conservation  of  energy  equation  becomes 

(H1  '  V  *  %  +  VS*  C‘  (B75) 

the  conservation  of  the  mass  of  water  vapor  is 

(MV,1  *  Vs* *  Vm.*0  •  (876) 

and  the  conservation  of  dry  air  mass  Le  . 


Introducing  the  quantities  AQ,*  M^q  o 

AQ  -  Hx  +  Qm  +  Qt  -  Q0 

Vo  =  Vl  +  Vm 


Hg  AQ 


(B77) 

(B78) 

(B79) 

(B80) 

(B8l) 


V,2  *  Vo  v.  ,  .  •  .  >  (B82) 

•\,2  =  Vo  -  -  (B83) 

All  of  the  other  relationships  remain  the  same  as  in  the  derivation  of  the 
transient  analysis..  The  shelter's  dry-bulb  temperature  is  still  define  by 
Equation  (B20)  and,  its  relative  humidity  is  still  defined  by  Equation  (plf-i)  with 
the  exception  that  the  terms  AQ,  My  q>  and  Ma^0  must  be  defined  as  in  Equations 
(B78)  through  (b8o) .  One  consequence  of  the  steady-state  derivation  is 
shelter  deminsions  and  volume  do  not  enter  into  the  calculations  of  "the  shelter ' s 
psychrometric  condition. 
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